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Abstract. Models of a Cen A & B have been computed using the masses determined by Pourbaix et al. (20021 
and the data derived from the spectroscopic analysis of Neuforge and Magain 1199711 . The seismological data 
obtained by Bouchy and Carrier (2001 2002) do help improve our knowledge of the evolutionary status of the 
system. All the constraints are satisfied with a model which gives an age of about 6 Gyr for the binary. 
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1. Introduction 

The binary system Alpha Centauri offers a unique oppor- 
tunity to test our knowledge of stellar physics in solar- type 
stars other than the Sun. On the one hand, its proximity 
allows very good determinations not only of its parallax 
but also of its orbital parameters so that the masses and 
the luminosities of both components can be determined. 
On the other hand, effective temperatures and chemical 
compositions can be obtained through spectroscopic anal- 
yses. 

A number of theoretical calibrations of this system 
have already been carried out. The reader will find ref- 
erences to earlier works in the publications of Guenther 
and Demarque (2000) and Morel et al. (2000). Very re- 
cently, a new and very precise analysis by Pourbaix et 
al. (2002) has led to new and strong constraints on the 
masses. Alpha Cen A thus becomes an excellent candidate 
for asteroseismological studies since solar-like p-mode os- 
cillations have recently been discovered from ground-based 
analysis (Bouchy and Carrier 2001, Carrier et al. 2002, 
Bouchy and Carrier 2002). 

In a very recent calibration and seismic analysis, 
Thevenin et al. (2002) have found that the asteroseismo- 
logical constraints, namely the large and small frequency 
spacings, required a slight decrease in the mass assigned 
by Pourbaix et al. (2002) to a Cen A. Their model repro- 
duces the large and small frequency spacings although the 
fit in the frequencies themselves is not achieved. 

Using an entirely new evolution code,which does not 
yet include diffusion of chemical elements, we find models 
for the a Cen binary system which agree not only with 
all the non-asteroseismic constraints, including the masses 



and the requirement of the same age for both components, 
but also with the astcroscismic constraints. 



2. Observational constraints 
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Using data from different sources, Pourbaix et al. (2002) 
have recently improved the precision of the orbital param- 
eters of a Cen. They have adopted the parallax derived 
by Soderhjelm ([1999), vo = 747.1 ± 1.2 mas, and give for 
the masses M A /M Q = 1.105 ± 0.0070 and M B /M Q = 
0.934 ±0.0061. 

From their spectroscopic analysis of the a Centauri 
system, Neuforge and Magain (1997) have determined the 
effective temperature, the gravity and the metallicity of 
both components. The effective temperatures are given 
by T cS:A = 5830 ± 30K and T e ff,i? = 5255 ± 50K. The 
gravities are given by loggA — 4.34 ± 0.05 and loggs = 
4.51 ±0.08. Since the luminosity can be expressed in terms 
of the effective temperature and of the gravity through the 
relation 

log l/Lq = log m/m q + 4io g r cff /r effi0 - io gg /g Q , 

the error box in the (logT e fT, log g) plane appears as a 
parallelogram in the HR diagram. 

For the metallicity, they obtained [Fe/H]^ = 0.25 ± 
0.02 and [Fe/H] B = 0.24 ± 0.03. To deduce Z/X, we as- 
sume that this ratio is proportional to the abundance ra- 
tio Fe/H and we adopt the solar value (Z/X)q — 0.023 
with an uncertainty of 10 %, given by Grevesse and 
Sauval (1998}. We then have the same value for both 
components of a Cen: Z/X = 0.040 ± 0.005. However, 
Grevesse and Sauval l|2002|) now favour a lower value 
(Z/X)q — 0.0209 with the same uncertainty. With this 
solar value, we would have Z/X = 0.037 ±0.004. It seems 
safe to say that Z/X is between 0.033 and 0.045. 
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The above non-asteroseismic constraints are summa- 
rized in Tabled 

o.c 

Table 1. Non-asteroseismic constraints 





A 


B 


M/Mq 


1.105 ±0.0070 


0.934 ± 0.0061 


log T cff 


3.7657 ± 0.0022 


3.7206 ± 0.0041 


logs 


4.34 ± 0.05 


4.51 ± 0.08 


z/x 


0.039 ± 0.006 



3. The models 

A number of evolutionary sequences have been com- 
puted from the main sequence with CLES (Code Liegeois 
d'Evolution Stellaire). We use the CEFF equation of state 
(Christensen-Dalsgaard and DaDDen ri992[> . The opacities 
are those of the Lawrence Livermore National Laboratory 
(Iglesias and Rogers I1996fl completed with the Alexander 
and Ferguson l|1994f) opacities at low temperature, both 
tables being smoothly joined in the temperature range 
defined by 3.95 < logT < 4.10. The nuclear reaction rates 
are from Caughlan and Fowler l|1988[) and the screening 
factor comes from Salpeter l|1954[) . The boundary condi- 
tions at the photosphere are deduced from Kurucz £1998) 
models. We use the mixing length theory of convection. 
The convection parameter a = i/Hp is left free and is 
used to force the evolutionary track through the error box 
around the current position of the star in the HR diagram. 
The suitable value of a depends on the initial chemical 
composition. It turns out that, for a given composition, 
we can always use the same value of a for both compo- 
nents, so that the equality a a = cub comes out naturally 
and is not to be considered as a supplementary constraint. 
The code does not include diffusion. 

We have computed a number of models for both com- 
ponents, with different initial chemical compositions, de- 
fined by (X, Z), and convection parameter a. We require 
that both components of the binary system reach their 
respective positions in the HR-diagram at the same age. 
This requirement of simultaneity determines a line in the 
(X, Z)-diagram to the right of which the B component 
reaches its observed position in the HR-diagram after the 
A component has already left its own observed position 
(see Fig. QJ. In addition, the chemical evolution of the 
Galaxy suggests a correlation between X and Z. The re- 
lation X = 0.76 — 3Z which has been used by Pols et 
al. (1998) is also shown in Fig.^ Although such a relation 
may have a statistical significance, it may be hazardous to 
strictly adhere to it in a particular case. These constraints, 
together with the constraint on the Z/X value, delimit a 
permitted area in the (X, Z) diagram. 

4. Asteroseismology of a-Centauri 

Bouchy and Carrier l|2001|) . Carrier et al. (|2002|) . and 
Bouchy and Carrier j2DEl) have recently observed the os- 



no simultaneity 
in error boxes 



Fig. 1. Constraints on the initial chemical composition 
coming from the spectroscopy (solid line) and from the 
simultaneity condition (dashed line). The dotted line rep- 
resents the chemical evolution of the Galaxy. The crosses 
indicate the different compositions for which evolution se- 
quences have been computed. The black square represents 
the model discussed in detail in this paper. 

Table 2. Mode frequencies (in /iHz) of a Cen A. 





Observations" 




Our model 






I = 


I = 1 


I = 2 


I = 


I = 1 


I = 2 


n — 


15 






1833.1 


1730.5 


1778.5 


1828.0 


n — 


16 


1841.3 


1887.4 


1934.9 


1834.9 


1882.8 


1932.6 


n — 


17 




1991.7 


2041.5 


1939.3 


1987.7 


2038.1 


n — 


18 




2095.6 


2146.0 


2044.4 


2093.7 


2144.7 


n — 


19 


2152.9 


2202.8 


2251.4 


2150.8 


2200.4 


2251.7 


n — 


20 


2258.4 


2309.1 


2358.4 


2257.5 


2307.5 


2358.7 


n — 


21 


2364.2 


2414.3 


2464.1 


2364.2 


2414.4 


2465.9 


n — 


22 


2470.0 


2519.3 


2568.5 


2471.0 


2521.5 


2573.1 


n — 


23 


2573.1 


2625.6 




2578.0 


2629.0 


2681.0 


n — 


24 


2679.8 


2733.2 


2782.9 


2685.5 


2736.9 


2789.2 


n — 


25 


2786.2 


2837.6 


2887.7 


2793.4 


2845.2 


2897.6 



a Bouchy and Carrier (|![002). 

dilation spectrum of a Cen A with the CORALIE spec- 
trograph. In the first paper, they observed the star during 
5 nights, and were able to identify some / = and I = 1 
modes, but the resolution was too low to identify sepa- 
rately the / = and the I = 2 modes. In their next run 
of observations, they observed the star during 11 nights, 
and were able to identify many I = 0, I = 1 and I = 2 
modes. They could not resolve the I = 3 modes. Their 
latest results are summarized in Table FJ1 They obtain an 
uncertainty of 0.46/iHz on these frequencies, but add that 
an error of ±1.3/iHz could have been introduced for some 
of the frequencies due to aliases. For this reason, we will 
adopt this last value for the observational error bars. 

We have calculated the oscillation frequencies of 
a Cen A using a standard adiabatic code (Boury and 
al.^975), for a grid of models where the mass Ma is fixed 
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at the value determined by Pourbaix et al. J2002J), i.e., 
Ma = 1.105Mq. The only free parameters are a, X, and 
Z, and the last two can only vary within the permitted 
area illustrated in Fig.^ The models we have explored are 
shown as crosses in Fig. ^ The models which fall outside 
of the permitted area were explored for reason of com- 
pleteness only. 

We calculate the modes I = 0, I = 1 and I — 2 for n = 
15 to 25. We then calculate the large frequency spacings 
&-vi,n = W,n — vi,n-x and the small frequency spacings 
Svo,n = W,n — vi + 2,n~i, and their averages over n = 15 to 
25.' 

5. Discussion and Conclusions 





3.78 3.77 



3.74 

'°0 T 8lf 



Fig. 2. Evolutionary tracks in the HR diagram for the 
sequences obtained with X — 0.70, Z — 0.0275, and a = 
1.8. The 1 a error boxes on log L/Lq and logT c ff for both 
components appear as parallelograms. The crosses denote 
the model retained for our discussion. 



Table 3. Our model 



A 



M/M & 


1.105 0.934 


L/L Q 


1.54 0.544 


T c b 


5821K 5291K 


X 


0.70 


z 


0.0275 


a 


1.8 


age 


6.41Gyr 



We will discuss the results corresponding to one of our 
best solutions (see Table shown as a black square in 
Fig. n Let us notice that for the same physics, a solar 
calibration gives a = 1.77. The trajectories in the HR 
diagram for the two components of the system are shown 
in Fig.0 The cross on the A component track in the HR- 
diagram shows the model retained for our discussion. We 
list the frequencies we have obtained for this particular 




v n | (mHz) 

Fig. 3. Large frequency spacing as a function of the fre- 
quency. The symbols indicate the values obtained by 
Bouchy and Carrier i|2002[) . with error bars of ±1.3/iHz. 
The three lines correspond to the model shown as a cross 
in Fig.©, for I = 0, I = 1 and I = 2. 




v n | (mHz) 

Fig. 4. Small frequency spacing 5v$ as a function of the 
frequency. The symbols represent the values obtained by 
Bouchy and Carrier H2()02|l . with error bars of ±1.3/iHz. 
The line corresponds to the model shown as a cross in 
Fig.0 

model in Table [2 Finally, we show in Figures and 0] the 
results for the large and the small frequency spacings. The 
age of this particular model is t = 6.4 Gyr. The position 
of the B component in the HR diagram at the same age 
is also marked with a cross in Fig. [2J 

This model satisfies all the non-asteroseismic con- 
straints listed in Tabled Moreover, the two components of 
the system have the same age. As to the asteroseismic con- 
straints, the average of the large spacings is slightly larger 
than the one obtained from the observations, but those 
large spacings remain within the error bars, as shown in 
Fig. |2| The values for the small spacings fall within the 
error bars as well, as shown in Fig.^J Finally, we note that 
this solution was found for a value of the mixing-length 
parameter a very close to the solar value. 
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The value of the average small spacing {6vq) is sensitive 
to the structure of the core of the star, and therefore to its 
chemical composition. Let us note that our retained model 
does not have a convective core, but it nevertheless leads 
to an excellent agreement for the small frequency spacings. 
The value of (Svo) increases as we lower the value of X. 
But all the models which lie on a curve more or less par- 
allel to the dashed curve in Fig. ^ provide the same value 
for (Svo). We are unable, given the resolution of the obser- 
vations available, to discriminate between those solutions. 
To go one step further we would need additional observa- 
tional data, such as the I = 3 modes, or even better, the 
oscillation frequencies of a Cen B. 

We have been able with our model to reproduce quite 
well, and certainly within the observational error bars, 
both the large and the small frequency spacings observed 
by Bouchy et al. with the CORALIE spectrograph for 
a Cen A, while respecting all the non-asteroseismic con- 
straints of the system. We must however be very careful 
with the interpretation of the results, since we have not 
included the diffusion of helium and other heavy elements 
in our stellar evolution code. It is well known that the 
inclusion of diffusion has greatly improved the agreement 
between solar models and observations (see, e.g., Bahcall 
and Pinsonneault 1995, Morel et al. 1997). However, sev- 
eral papers have examined the influence of diffusion on 
models of a Cen A (Edmonds et al. 1992, Turcotte and 
Christensen-Dalsgaard 1998) and they have shown that 
diffusion of helium has no detectable effect on the frequen- 
cies of the low ^-degree high-order p-modes, and therefore 
no effect on the large and the small frequency separations, 
in the frequency range of the spectrum of a Cen A. But 
it has a substantial effect on the age and on the initial 
composition. According to Edmonds et al. (1992) diffusion 
lowers the age of the star by about 10%, and it increases 
the initial helium abundance to Y=0.30. 

In their very recent paper, Thevenin et al. ( 2002) pro- 
pose a model for a Cen A which gives acceptable values for 
the large and the small frequency spacings, but they had 
to change the mass Ma of the star. In addition, their fre- 
quencies differ from the observed frequencies by « 30/iHz. 
It is well-known that the exact values of the frequencies 
themselves depend on the details of the star's atmosphere, 
where the pulsation is non-adiabatic. A discrepancy of 
30/iHz seems rather large, however, since typical discrep- 
ancies for solar p-mode frequencies are smaller than 10/iHz 
(see, e.g., Guenther et al. 1996, Guenther & Demarque 
1997). We show in Fig. the differences between the cal- 
culated frequencies and the observed frequencies for our 
model. We see that these differences display a pattern sim- 
ilar to what is found for the Sun (Guenther et al. 1996, 
Guenther & Demarque 1997). 

The evolution code used by Thevenin et al. (2002) 
differs from ours in two ways: they use the Canuto & 
Mazitelli (1991, 1992) theory of convection (hereafter CM) 
below an optical depth r = 20 while we use the basic 
Bohm-Vitense (1958) mixing-length theory of convection 
(hereafter MLT), and they include the effects of helium 
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Fig. 5. Differences between the calculated frequencies and 
the observed frequencies for our model 

diffusion. Morel et al. (2000) obtain very different ages 
for the a Cen system when using the CM or the MLT 
theories of convection: 4086 Myr and 2710 Myr respec- 
tively. In fact, changing the treatment of the convection, 
while keeping all other stellar parameters the same, will 
affect the value of the effective temperature at a given 
age (or luminosity). Therefore, in order to go through a 
given position in the HR-diagram, it is necessary to mod- 
ify, e.g., the helium abundance. Morel et al. (2000) had to 
use Y=0.284 in MLT models, and Y=0.271 in CM models. 
This of course changes considerably the age of the star. On 
the other hand, diffusion only changes the initial helium 
abundance by about 0.5% (Turcotte and Christensen- 
Dalsgaard 1998), and the age by 0.4Gyr (Edmonds et al. 
1992). Furthermore, the effects of microscopic diffusion 
could be even reduced if turbulent mixing is taken into 
account (Turcotte et al. 1998). In addition, Thevenin et 
al. (2002) use observations different from the ones we use 
to calibrate their models (Morel et al. 2000). The effective 
temperature they use for a Cen A is 40K lower than ours. 
This will also affect the luminosity. Moreover, they use 
a rather smaller value for the metallicity [Fe/H]=0.20 in- 
stead of our [Fe/H]=0.25. The discrepancies between the 
results of Thevenin et al. (2002) and the results we present 
here could result from those differences. It would of course 
be very useful to compare the results of different stellar 
evolution codes in a more consistent and systematic man- 
ner. 
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